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HIGHLIGHTS 


•  A  numerical  model  is  developed  for  a  H202-direct  borohydride  fuel  cell  (H202-DBFC). 

•  Global  reactions  in  the  model  are  calibrated  vs.  a  simple  experimental  H202-DBFC. 

•  The  calibrated  model  is  used  to  predict  the  losses  from  parasitic  side  reactions. 

•  Maximum  power  density  and  system  efficiency  coincide  due  to  parasitic  side  reactions. 

•  H+  reduction  contributes  to  cathode  current  density  when  H202  is  transport  limited. 
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A  numerical  model  with  global  reaction  rates  is  calibrated  to  measurements  from  a  simple  hydrogen- 
peroxide  direct-borohydride  fuel  cell  (H202-DBFC),  and  then  used  to  unravel  complex  electrochemical 
and  competing  parasitic  reactions.  In  this  H202-DBFC,  fuel  (1-50  mM  NaBH4/2  M  NaOH)  is  oxidized  at  a 
Au  anode  and  oxidizer  (10—40  mM  H202/1  M  H2S04)  is  reduced  at  a  Pd:Ir  cathode.  Polarization  curves 
and  electrode  potentials,  as  functions  of  fuel  and  oxidizer  feeds  support  global  reaction  rate  parameter 
fitting.  The  measurements  and  calibrated  model  showed  H202  decomposition  at  the  cathode  depresses 
open  circuit  voltage  from  3.01  V  theoretical  to  1.65  V,  and  when  H202  supply  is  limited,  cathode  po¬ 
tentials  are  sufficiently  negative  to  make  H+  reduction  to  H2  thermodynamically  favorable.  Calibrated 
model  results  show  that  thin  concentration  boundary  layers  limit  reactant  utilization  and  current  den¬ 
sity.  Decreasing  the  inlet  concentrations,  flow  rates,  and  cell  voltage  slow  parasitic  reactions  and  favor 
desirable  charge  transfer  reactions.  Peak  conversion  efficiency  and  peak  power  density  coincide  because 
thermodynamic  efficiency  and  parasitic  reaction  rates  decrease  (relative  to  charge  transfer  reaction  rates) 
with  increasing  current  density.  We  conclude  that  the  performance  of  a  fuel  cell  with  parasitic  side 
reactions  can  be  predicted  through  numerical  modeling. 

Published  by  Elsevier  B.V. 


1.  Introduction 

Direct  borohydride  fuel  cells  (DBFCs)  have  been  considered  for 
applications  ranging  from  portable  consumer  electronics  to  vehicle 
propulsion.  The  term  DBFC  refers  to  diverse  fuel  cell  configurations 
in  which  BH4  is  oxidized  at  the  anode;  review  papers  1-8]  sum¬ 
marize  the  variety  of  DBFC  chemistries,  catalysts  and  topologies 
explored  to  date.  Two  major  classes  include  those  utilizing  acidic  or 
alkaline  oxidizer  solutions.  The  DBFC  studied  herein  uses  NaBH4 / 
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NaOH  fuel  and  H2O2/H2SO4  oxidizer,  and  is  referred  to  herein  as  a 
H2O2-DBFC. 

This  H2O2-DBFC  offers  several  prospective  advantages.  First,  the 
reactants  can  be  stored  at  ambient  temperature  and  pressure  with 
concentrations  of  several  molar,  rather  than  as  a  high-pressure  gas 
or  cryogenic  liquid.  Second,  when  high  reactant  storage  concen¬ 
trations  are  combined  with  the  3.01  V  theoretical  cell  potential  for 
the  NaBH4/NaOH  -  H2O2/H2SO4  couple,  the  theoretical  energy 
density  and  specific  energy  of  the  reactants  surpass  those  of 
advanced  Li-ion  batteries  [2].  Finally,  DBFCs  that  consume  stored 
H2O2  oxidizer  can  operate  in  anaerobic  environments  [9,10]  (e.g. 
underwater  or  in  space)  and  offer  higher  power  density  than 
systems  that  rely  on  O2  due  to  the  facile  kinetics  of  H2O2 
reduction. 
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The  desired  electrochemical  reactions  occurring  in  a  H2O2-DBFC 
are  written  in  (R1 )— (R3).  The  net  cell  reaction  (R1  consumes  BH4 , 
H202,  OH-,  and  H+  and  forms  B02  and  H2O  for  a  standard  cell 
potential  =  3.01  V.  The  anode  and  cathode  half-cell  reactions 
are  in  (R2)  and  (R3),  respectively.  In  (R2),  BH4  and  OH-  are  elec- 
trochemically  oxidized  to  B02  and  water.  Electrons  provided  by 
(R2)  travel  through  an  external  circuit  to  be  consumed  by  the 
cathode  half-cell  electrochemical  reduction  of  H2O2  and  H+  (R3). 

BH4  +  80H-  +  4H202  +  8H+  BO2  +  14H20  (R1 ) 

BH4  +80H“^B0J  +6H20  +  8e-  E$2  =  -1.246  V  vs.  RHE 

(R2) 

H202  +  2H+  +2e-^2H20  E^3  =  1.763  V  vs.  RHE  (R3) 

Parasitic  chemical  reactions  can  also  occur;  they  uselessly 
consume  reactants  and  depress  cell  voltage  by  shifting  concentra¬ 
tions  near  the  electrodes.  At  the  anode  catalyst,  when  the  supply  of 
OH-  is  insufficient,  heterogeneous  hydrolysis  11]  converts  BH4  to 
B02  and  H2  via  parasitic  chemical  reaction  (R4). 

BH4  +  2H20  ->  BO2  +  4H2  (R4) 

Similarly,  if  the  supply  of  H+  in  the  oxidizer  solution  is  insuffi¬ 
cient,  H202  undergoes  heterogeneous  decomposition  at  the  cath¬ 
ode  catalyst  [12-14]  to  form  water  via  the  parasitic  chemical 

reaction  (R5). 

2H202— >2H20  +  02  (R5) 

Reactions  (R1)-(R5)  each  summarize  multiple  elementary  in¬ 
termediate  steps.  For  example,  transferring  8e-  to  the  anode  in 
reaction  (R2)  requires  at  least  8  elementary  reactions.  The  micro¬ 
kinetics  of  elementary  reactions  occurring  at  a  DBFC  anode  has 
been  explored  by  Rostimikia  et  al.  using  density  functional  theory 
[15].  The  rate  determining  elementary  reactions  and  related  kinetic 
parameters  at  the  anode  and  cathode  of  a  H202-DBFC  are,  however, 
still  uncertain. 

Experimental  investigations  of  H202-DBFCs  indicate  that  they 
do  not  meet  their  theoretical  promise  primarily  because  of  trans¬ 
port  limitations  in  the  liquid  reactant  flow  channels  and  the  loss  of 
fuel  and  oxidizer  via  parasitic  chemical  reactions  at  each  electrode 
[11,16-19].  Unfortunately,  the  interdependencies  of  the  H202-DBFC 
electrochemistry  and  parasitic  reactions  are  difficult  to  decouple 
experimentally  because  of  their  complexity.  Carefully  constructed 
numerical  models  are  well  suited  to  such  problems,  because  they 
can  capture  the  complexity  while  revealing  parameters  that  are 
inaccessible  experimentally.  We  previously  developed  a  numerical 
model  to  understand  the  relationships  between  cell  geometry, 
operating  conditions,  and  performance  in  an  ideal  H202-DBFC 
where  all  reaction  rates  were  fast  and  parasitic  reactions  were 
omitted  [20,21  ].  For  this  present  work,  we  revise  the  existing  model 
to  include  realistic  reaction  rates  and  parasitic  reactions.  Experi¬ 
ments  with  a  H202-DBFC  guide  selection  of  appropriate  global  re¬ 
actions  for  the  model,  and  provide  measurements  to  which  the 
model  is  calibrated.  The  model  is  calibrated  by  fitting  kinetic  terms 
in  the  global  reaction  rates  to  H202-DBFC  polarization  curves 
measured  over  a  range  of  BH4  concentrations.  The  simplified  re¬ 
action  mechanism  in  the  revised  model  reduces  the  large  number 
of  uncertain  kinetic  parameters  in  comparison  to  elementary 
microkinetic  models,  and  thereby  makes  fitting  to  measurements 
more  tractable.  With  this  simplification,  the  global  reaction  rates 
with  fitted  parameters  can  provide  sufficient  accuracy  for  cell-level 


studies  as  long  as  they  reasonably  capture  the  rate-determining 
processes. 

Following  calibration  against  the  experimental  H202-DBFC,  the 
model  is  used  to  investigate  the  effects  of  parasitic  reactions  on 
H202-DBFC  performance  by  examining  a  baseline  case  and  then 
departures  from  the  baseline  with  variations  in  fuel  inlet  concen¬ 
tration  and  flow  rate. 

2.  Model  development 

A  steady  state  isothermal  H202-DBFC  having  the  form  shown  in 
Fig.  1  was  modeled  with  a  2-D  finite  volume  approach.  Aqueous  fuel 
(1-50  mM  NaBH4/2  M  NaOH)  and  oxidizer  (10-40  mM  H202/1  M 
H2S04)  flow  through  rectangular  channels  separated  by  a  Nation 
117  membrane.  Catalyst-coated  channel  walls  constrain  the  flows 
and  act  as  electrodes.  BH4  is  oxidized  at  the  Au  anode  and  H202  is 
reduced  at  the  Pd:Ir  cathode.  Na+  cations  cross  the  membrane  to 
charge  balance  the  cell,  bringing  H20  with  them  via  electro-osmotic 
drag.  The  geometry  in  Fig.  1  represents  a  class  of  DBFCs  having 
separated  electrodes  and  membranes.  Such  DBFCs  have  been 
considered  in  earlier  studies  [22,23]  because  they  are  simple  to 
fabricate,  resist  precipitate  accumulation,  expose  the  full  mem¬ 
brane  to  the  aqueous  electrolytes,  and  orient  the  electric  field  so 
that  it  aids  reactant  transport  [21  .  This  DBFC  configuration  can  be 
represented  in  a  computationally  efficient  2D  model,  when  the 
channel  sidewalls  are  inert  (electrochemically  inactive)  and  spaced 
widely  enough  to  have  negligible  effect  on  the  hydrodynamics  of 
the  reactant  flows.  The  model  predicts  cell  current  when  the  cell 
geometry,  inlet  flow  rates,  inlet  concentrations,  and  cell  voltage  are 
specified. 

Only  a  summary  of  the  transport  calculations  and  solution 
approach  are  provided  herein  because  they  have  been  described 
previously  [20,21  .  The  reaction  mechanisms  and  rates  are  dis¬ 
cussed  in  detail  because  they  are  specific  to  this  study. 

2  A.  Transport  in  the  channels  and  membrane 

State  variables  describe  the  reactant  solutions  for  computational 
cells  in  each  flow  channel:  pressure  P,  velocities  vx  and  vy,  electric 
potential  c[>,  and  species  mass  fractions  Y/<.  For  the  fuel  solution 
species,  k  =  BH4,  B02 ,  Na+,  OH-,  H2,  or  H20.  For  the  oxidizer  so¬ 
lution,  k  =  H202,  H+,  Na+,  SO42,  H2,  02  or  H20.  The  electrode¬ 
channel  interfaces  have  state  variables  0  and  Y/<,  and  the  mem¬ 
brane-channel  interfaces  have  state  variables  0,  Y/<  and  P.  The  value 
of  each  channel  state  variable  is  found  by  enforcing  the  associated 
conservation  equation;  conservation  of  mass  (P),  conservation  of 
species  (Y/<),  conservation  of  momentum  (v)  and  electroneutrality 
(0).  Solving  for  0  by  enforcing  electroneutrality  implies  that  elec¬ 
trochemical  double  layers  (the  only  regions  in  which  substantial 
deviation  from  electroneutrality  occur)  are  thin  and  effectively  part 
of  the  electrode-channel  interfaces.  The  interface  values  of  P,  0,  and 
Y/<  are  found  by  mass,  charge  and  species  flux  matches,  respectively. 
Fluxes  due  to  electrode  reactions  and  channel  transport  must 
match  at  electrode-channel  interfaces,  and  fluxes  due  to  membrane 
transport  and  channel  transport  must  match  at  the  mem¬ 
brane-channel  interfaces.  Mass,  mole  and  charge  densities  are 
estimated  using  apparent  molar  solute  volumes  in  the  reactant 
solutions  [21  . 

Boundary  conditions  at  the  inlets  and  outlets  include  fully 
developed  momentum  boundary  layers  (9vx/0x  =  0,  dVy/dy  =  0)  and 
no  electric  field  in  the  x-direction  (00/9x  =  0).  Mass  fractions  are 
specified  at  the  inlet  (dictated  by  fuel  and  oxidizer  concentrations) 
and  pressure  is  specified  at  the  outlet  (P  =  0).  Boundary  conditions 
at  the  electrodes  include  no-slip  and  impermeable  boundary 
(vV  =  0,  a P/ay  =  0).  The  anode  is  assigned  0  =  0,  making  it  the 
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Fig.  1.  Illustration  of  our  DBFC  cell  configuration  showing  the  cell  geometry  and  major  species  flows,  reactants,  and  products. 


electric  potential  reference  for  the  cell,  and  the  cathode  is  assigned 
f  =  Vceii.  The  only  boundary  condition  at  the  membrane  is  no-slip 

(avx/a*  =  o). 

Several  simplifications  are  appropriate  given  the  cell  geometry 
and  operating  conditions.  The  reactant  flows  are  laminar 
(. Re  «  2300),  the  cell  is  isothermal  (heat  capacities  of  reactant 
solutions  are  small  compared  to  heating  rate),  the  membrane  is 
fully  hydrated  (in  contact  with  aqueous  solutions  on  both  sides) 
and  there  are  no  homogenous  reactions  in  the  channel.  All  re¬ 
actions  occur  at  the  catalyst/fluid  interface.  The  reactant  solutions 
are  assumed  to  be  ideal  (i.e.  activity  coefficients  yk  =  1 )  because  the 
activity  coefficients  for  reactant  species  not  in  excess  (BH4  and 
H2O2)  vary  little  with  concentration  over  the  ranges  examined  in 
this  study  [24  .  When  dyk/dCk  ~  0,  non-idealities  have  little  influ¬ 
ence  on  the  rates  of  transport  [25]. 


2.2.  Electrode  reactions 

The  net  rate  of  each  charge  transfer  reaction  at  an  electrode  is 
modeled  as  the  sum  of  anodic  and  cathodic  rates  as  in  Eq.  (1), 
where  the  anodic  direction  provides  electrons  to  the  electrode 
and  the  cathodic  direction  removes  them.  The  net  rate  depends 
on  rate  constants  k ,  concentrations  Ck,  stoichiometric  coefficients 
v,  electric  potential  difference  across  the  electrode-channel 
interface  A0,  symmetry  factors  fi,  number  of  electrons  trans¬ 
ferred  in  the  rate  limiting  step  ne,  and  /  =  F/RT.  Net  rates  for 
chemical  reactions  omit  the  electric  potential  activation  energy 
barriers  by  setting  ne  =  0;  in  this  case  the  anodic  and  cathodic 
rates  simplify  to  forward  and  reverse  rates  as  in  Eq.  (2).  The 
nomenclature  for  symbols  appearing  in  equations  is  summarized 
in  Appendix  A. 


Eq.  (3)  adjusts  the  reaction  fluxes  based  on  the  ratio  of  the  active 
electrochemical  surface  area  to  the  geometric  electrode  area. 

J  k, rxn  =  ®  'y  ^  rqvk,q  (3) 

q 

Calibration  experiments  are  used  to  fit  rate  parameters  ka  and  pa 
for  each  electrochemical  reaction,  /</for  each  chemical  reaction,  and 
£  for  each  electrode.  Cathodic  rate  constants  kc  are  chosen  to 
maintain  thermodynamic  consistency  given  the  fitted  values  of  ka. 
Thermodynamic  consistency  in  this  case  is  defined  as  zero  net  rate 
for  each  reaction  under  standard  conditions  (1  M  concentration, 
298  K  temperature,  and  101.3  kPa  pressure)  at  its  equilibrium 
electrode  potential.  This  is  equivalent  to  solving  kc  =  kae^° . 

To  find  values  of  ka,  (3a  and  £  that  agree  with  the  calibration 
experiments,  a  comparison  function  is  developed.  The  comparison 
function  calls  the  H2O2-DBFC  model  code  to  generate  a  predicted 
polarization  curve,  and  then  calculates  the  difference  (residual) 
between  measured  and  predicted  polarization  curves  at  each  cell 
potential.  The  MATLAB  function  Isqnonlin  with  a  Newton  search 
algorithm  minimizes  the  L2-norm  of  the  residuals  by  adjusting  the 
fit  parameters,  ka,  j3a,  and  £.  The  fit  parameters  are  constrained  to 
intervals  (0,oo),  (0,1]  and  [1,30]  respectively.  Trial  and  error  with  fit 
parameter  initial  guesses  show  that  the  initial  guess  must  produce 
a  polarization  curve  differing  from  the  measurements  by  less  than 
one  order  of  magnitude  or  the  fitting  algorithm  would  likely  find  a 
local  minimum.  At  least  three  widely  spaced  starting  guesses  are 
evaluated  in  each  fitting  effort,  and  the  fit  is  not  accepted  as  “final” 
until  all  three  guesses  result  in  the  same  fit,  suggesting  the  fit  may 
be  global. 

2.3.  Calculating  H2O2-DBFC  performance  metrics 


r  —  ka  n  C  eUefja^<p 

k 


ne(3j  A0 


r 


Cathodic  symmetry  factors  (3C  are  estimated  as  /3C  =  1  ~/3a  and  ne 
is  assumed  to  be  1  based  on  the  expectation  that  one  elementary 
charge-transfer  step  (with  single  electron  transfer)  is  likely  rate- 
limiting  for  the  overall  reaction  [26,27  .  Concentrations  Ck  near 
the  electrode  and  electric  potential  differences  A  f  are  provided  by 
the  transport  portion  of  the  model.  __ 

The  net  flux  of  each  species  k  from  an  electrode,  J  krxn,  is  the 
sum  of  contributions  from  all  reactions  q  at  the  electrode  which 
involve  species  k,  as  shown  in  Eq.  (3).  A  roughness  parameter  £  in 


The  H2O2-DBFC  performance  is  characterized  herein  by 
several  metrics:  reactant  utilization,  coulombic  efficiency,  ther¬ 
modynamic  efficiency,  coulombic  utilization,  and  conversion  ef¬ 
ficiency.  These  metrics  describe  the  effectiveness  with  which 
reactants  are  used  on  a  single  pass  through  the  cell,  and  the 
prospective  performance  improvement  provided  by  recircula¬ 
tion,  making  them  relevant  to  both  single-pass  and  recirculated 
H2O2-DBFCS.  Furthermore,  these  metrics  are  strongly  influenced 
by  parasitic  side  reactions,  which  are  present  in  both  single-pass 
and  recirculated  H2O2-DBFCS. 

Reactant  utilization  7]rUtk  is  the  fraction  of  reactant  k  consumed 
before  the  bulk  flow  carries  it  out  of  the  cell.  It  is  defined  for  the  cell 
geometry  in  this  study  by  Eq.  (4),  where  Lceii  is  the  channel  length 
from  inlet  to  outlet,  J  k  elec  is  the  local  flux  of  reactant  k  to  the 
electrode,  Ck\ n  is  the  concentration  of  k  at  the  channel  inlet  and  V[n 
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is  the  solution  volumetric  flow  rate  at  the  channel  inlet.  Fig.  1  shows 
that  the  x-coordinate  runs  from  inlet  to  outlet. 


Vm  ,k 


Coulombic  efficiency  pce<k  is  the  fraction  of  consumed  reactant  k 
which  contributes  to  the  desired  charge  transfer  reactions,  and  is  a 
measure  of  losses  to  parasitic  side  reactions.  It  is  calculated  via  Eq. 
(5),  where  i  is  local  current  density  and  ve  and  vk  are  the  stoichio¬ 
metric  coefficients  for  electrons  and  species  k  in  the  charge  transfer 
reaction(s),  respectively.  Average  coulombic  efficiency  over  the 
whole  channel  is  fjcek. 


Vce 


i 


J  k,e\ecve 


The  thermodynamic  efficiency,  rjth.  is  the  ratio  of  useful  (elec¬ 
trical)  work  done  by  the  fuel  cell  to  the  free  energy  change  of  the 
reactants.  Thermodynamic  efficiency  is  less  than  1  due  to  activation 
losses  at  each  electrode,  mass  transport  and  ohmic  losses  in  each 
channel,  and  ohmic  losses  in  the  membrane.  These  losses  are 
manifested  as  departures  from  the  theoretical  cell  voltage  (over¬ 
potentials),  so  the  thermodynamic  efficiency  can  be  written  as  in 
Eq.  (6). 


The  reactant  utilization,  coulombic  efficiency,  and  thermody¬ 
namic  efficiency  can  be  combined  to  form  two  other  performance 
metrics.  Coulombic  utilization  {vCuj<  =  Vru,kVceJ<)  is  the  fraction  of 
reactant  k  flowing  into  the  cell  which  contributes  to  charge  transfer 
reactions.  Conversion  efficiency  for  reactant  k  (r]c  k  =  VthVce.k)  is  the 
fraction  of  the  chemical  energy  carried  by  reactant  k  into  the  H2O2- 
DBFC  that  is  converted  to  electrical  energy.  Conversion  efficiency, 
which  depends  on  cell  operating  conditions  and  geometry,  de¬ 
scribes  the  amount  of  useful  electrical  energy  that  can  be  produced 
by  a  given  quantity  of  reactant  on  a  single  pass  through  the  H2O2- 
DBFC. 


3.  Experimental  procedures 

3  A.  H2O2-DBFC  fabrication 

A  H202-DBFC  matching  the  geometry  of  Fig.  1  was  designed  and 
fabricated.  The  cell  consisted  of  two  flat  graphite  electrodes  sepa¬ 
rated  by  flow  channel  masks  and  a  Nation  117  membrane  (Fig.  2). 
The  flow  channel  masks  were  0.5  mm  thick  PTFE  with  5  mm  wide 
channels  cut  in  them.  The  channels  revealed  a  reactant  port  at  each 
end,  so  that  when  the  cell  was  assembled,  reactants  flowed  through 
the  channels.  The  regions  of  the  graphite  plates  exposed  to  the  flow 
channels  each  had  an  electrodeposited  catalyst  layer.  Each  graphite 
electrode  was  in  contact  with  a  gold-coated  brass  current  collector. 
Holes  through  the  entire  assembly  admitted  bolts  to  compress  and 
seal  the  cell;  low  reactant  pressures  simplified  sealing.  The 
compliance  and  hydrophobicity  of  the  PTFE  channel  masks  and  the 
smoothness  of  the  graphite  plates  together  discouraged  fluid  entry 
between  the  layers.  Large  channel  mask  areas  ensured  dimensional 
stability  in  the  direction  perpendicular  to  the  graphite  plates  so  that 
the  channel  depth  was  dictated  by  the  mask  thickness.  Channels 
were  as  wide  as  possible  while  supporting  the  membrane  in  the 
gap  between  the  graphite  plates.  Components  were  selected  to 
resist  chemical  attack  by  the  corrosive  reactant  solutions.  Aside 
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Fig.  2.  Solid  model  of  the  DBFC  used  for  calibration  experiments.  Fuel  and  oxidizer 
flow  along  the  catalyst  strip. 


from  the  Nation  membrane  and  graphite  electrodes,  all  wetted 
components  were  made  of  PTFE  or  PVDF. 

Graphite  electrodes  were  prepared  for  catalyst  electro¬ 
deposition  by  rinsing  with  methanol  and  acetone,  drying,  and 
secondary  rinsing  with  18  MQ  water.  The  electrodes  were  sanded 
with  coarse  (P80  grit)  AI2O3  sandpaper  and  then  fine  (P600  grit)  SiC 
sandpaper,  and  resulting  graphite  powder  was  removed  by  soni- 
cation  in  an  Ultrasonic  Power  Corp.  model  5300/50-26-459  soni- 
cation  bath  for  20  min  at  200  power  (in  bags  with  18  MQ  water). 
Pre-  and  post-sonication  micrographs  showed  that  graphite  dust 
was  thoroughly  removed. 

The  electrodeposition  process  was  adapted  from  those 
described  by  Miley  28],  de  Leon  [29  ,  Urian  [22]  and  Bessette  [30]. 
Anode  and  cathode  catalysts  were  deposited  separately;  in  each 
case  the  cell  was  assembled  with  the  electrode  to  receive  the 
catalyst  (working  electrode,  or  WE),  a  graphite  counter  electrode 
(CE),  a  channel  mask,  and  no  membrane.  The  anode  deposition 
solution  consisted  of  2  mM  NaAuCU  and  200  mM  NaCl  in  18  MQ 
water.  The  cathode  solution  consisted  of  2  mM  PdCl2,  2  mM 
Na3lrCl6,  200  mM  KCl  and  0.1  M  HC1  in  18  MQ  water.  Each  was 
prepared  using  Alfa-Aesar  Premion  reagent  grade  salts.  Solutions 
were  circulated  through  the  cell  by  a  peristaltic  pump.  An  Autolab 
PGSTAT30  potentiostat/galvanostat  repeatedly  drove  current  pul¬ 
ses  (5  mA  cm-2  for  2.5  s  and  0  mA  cm-2  for  2  s)  through  the  cell 
while  monitoring  the  WE  potential  with  respect  to  a  Ag/AgCl 
reference  electrode  (RE)  in  the  beaker.  The  solution  residence  time 
in  the  cell  was  ~1  s,  so  this  procedure  replaced  the  solution  between 
pulses.  The  number  of  pulses  was  chosen  to  give  ~20  mg  cm~2 
coverage  (assuming  1:1  ratio  for  the  Pd  and  Ir  deposition  rates).  The 
inlet  and  outlet  ports  on  the  cell  were  reversed  half-way  through 
deposition  to  mitigate  layer  thickness  variation  from  inlet  to  outlet. 
Deposition  took  place  at  room  temperature  (~23  °C).  The  pro¬ 
cedures  and  results  for  electrochemical,  chemical,  and  morpho¬ 
logical  characterization  of  the  electrodes  are  provided  as 
Supplemental  information. 

The  Nation  117  membrane  was  prepared  by  boiling  in  3%  H2O2 
for  1  h  to  remove  surface  contaminants,  then  boiling  in  1%  NaOH  for 
1  h  to  convert  to  the  Na+  form,  and  finally  rinsing  with  18  MQ 
water. 

3.2.  H2O2-DBFC  test  stand  setup 

A  test  stand  was  built  to  control  the  flow  of  reactants  through 
the  cell  and  acquire  data.  The  reactants  were  prepared  as  described 
in  the  Supplemental  information  section,  and  reactant  solution 
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mass  densities  were  measured  to  evaluate  the  accuracy  of  model 
predicted  mass  densities.  The  pH  of  reactant  and  effluent  solutions 
was  measured  to  determine  whether  H+  and  OH-  crossover  must 
be  included  in  the  membrane  model.  Reactants  began  in  carboys, 
and  then  flowed  through  the  cell  and  an  adjustable  constriction 
used  to  control  the  flow  rate.  A  combination  of  gravity  and  17  kPa 
N2  overpressure  was  used  to  drive  the  flows,  ensuring  consistent 
flow  rate.  The  cell  was  operated  with  once-through  (not  recircu¬ 
lated)  flows  so  that  the  inlet  concentrations  were  known.  Further 
details  are  given  in  the  Supplemental  information  section. 

A  MACCOR  Model  2300  test  system  was  used  to  control  and 
measure  the  cell  voltage  and  current,  and  to  record  each  electrode 
potential  with  respect  to  a  Ag/AgCl  reference  electrode.  MACCOR 
electric  potential  measurements  were  calibrated  against  an  Autolab 
PGSTAT30  prior  to  the  experiments.  Cell  temperature  was 
measured  by  a  K-type  thermocouple  affixed  to  the  anode  current 
collector.  The  same  fuel  and  oxidizer  flow  rates  were  used  for  all 
experiments;  they  were  measured  twice  per  experiment  by 
measuring  the  time  to  fill  250  mL  graduated  cylinders.  Compiling 
all  flow  rate  measurements  gave  10.17  ±  0.46  mL  min-1  for  the  fuel 
and  10.05  ±  0.16  mL  min1  for  the  oxidizer.  The  time  measurement 
error  was  insignificant  compared  to  the  total  time  of  ~25  min,  so 
flow  rate  uncertainty  depended  on  the  graduated  cylinder  accuracy 
and  flow  rate  variation. 

3.3.  Electrochemical  evaluation  of  the  H2O2-DBFC 

Polarization  curves  were  measured  by  permitting  the  cell  to 
equilibrate  at  open  circuit  and  then  stepping  the  cell  potential  from 
0.4  V  up  to  open  circuit  while  measuring  the  current.  Consistent 
results  were  obtained  when  the  first  step  was  held  for  5  min  to 
thoroughly  reduce  the  cathode.  Each  subsequent  cell  potential  was 
held  for  2.5  min.  The  cell  rapidly  (<0.5  min)  reached  steady  state 
after  each  cell  potential  change.  Data  were  recorded  at  1  Hz  and  an 
average  current  at  each  voltage  step  was  compiled  by  averaging  the 
last  2  min  of  the  step.  Three  polarization  curves  were  measured  for 
each  test  condition  and  later  compiled  into  an  average  curve  with 
error  bars  marking  one  standard  deviation. 

Each  set  of  polarization  curve  measurements  began  and  ended 
with  a  baseline  measurement  to  reveal  changes  in  the  cell  state  (for 
example,  catalyst  oxidation  state).  Baseline  inlet  concentrations 
were  10  mM  NaBH4/2  M  NaOH  fuel  and  40  mM  H202/1  M  H2S04 
oxidizer.  This  combination  with  1:1  stoichiometry  (assuming  the 
ideal  reaction  (Rl))  was  chosen  to  make  the  baseline  curves  sen¬ 
sitive  to  changes  at  either  electrode.  Baseline  curves  varied  little, 
indicating  consistent  cell  state.  Solution  pH  measurements  before 
and  after  flowing  through  the  H202-DBFC,  and  solution  density 
measurements  are  provided  as  Supplemental  information. 

4.  Results  and  discussion 

4.1.  H2O2-DBFC  experiments 

The  five  polarization  curves  in  Fig.  3  (symbols)  are  measured 
with  the  H202-DBFC  operating  at  1,  2.5,  5, 10,  and  20  mM  BH4/2  M 
NaOH  fuel  and  40  mM  H202/1  M  H2S04  oxidizer.  The  variant  BH4 
concentrations  span  the  transition  from  fuel  limited  to  oxidizer 
limited  stoichiometry  with  respect  to  (Rl).  The  fuel-limited  case  is 
studied  with  1,  2.5  and  5  mM  BH4;  equivalent  fuekoxidizer  at 
10  mM  BH4;  and  an  oxidizer  limited  case  with  20  mM  BH4.  The 
reactant  concentrations  are  lower  than  those  in  previously  reported 
H202-DBFC  experiments  [22,31,32]  to  minimize  gas  formation  from 
BH4  hydrolysis  and  H202  decomposition  and  thereby  maintain  the 
validity  of  this  single-phase  flow  model.  The  lower  concentrations 
lead  to  lower  power  densities  than  in  prior  experiments.  A 


Fig.  3.  Polarization  curves  with  varying  BH4  concentration/2  M  NaOH  fuel,  40  mM 
H202/1  M  H2S04  oxidizer,  and  both  flow  rates  at  10  mL  min”1.  Symbols  are  measured 
and  curves  are  predicted  after  calibrating  rate  parameters  to  the  measurements.  Solid 
curves  exclude  reaction  R6  and  dash  curves  include  reaction  R6. 

qualitative  assessment  of  experimental  data  in  Fig.  3  indicates  that 
BH4  transport  limits  appear  in  the  fuel-limited  1, 2.5  and  5  mM  BH4 
curves  with  the  current  density  tapering  off  at  high  current  density, 
but  curves  measured  in  the  oxidizer  limited  cases  with  10  and 
20  mM  BH4  do  not  show  the  sharp  drops  in  current  density  ex¬ 
pected  for  H202  transport-limited  behavior.  This  lack  of  H202 
transport-limited  behavior  in  the  polarization  curves  is  observed 
despite  the  greater  diffusivity  of  BH4  (2.42  x  10-9  m2  s_1) 
compared  to  H202  (1.49  x  10-9  m2  s_1),  which  should  make  the 
H202  transport  limit  apparent  at  lower  BH4  concentration. 

The  lack  of  an  apparent  H202  transport  limit  in  Fig.  3  is 
explained  by  H+  reduction  occurring  at  the  cathode,  which  con¬ 
tributes  to  the  total  current  density  and  has  an  onset  potential  that 
causes  it  to  overlap  with  and  obscure  the  H202  transport  limit.  H+ 
reduction,  written  as  reaction  (R6),  is  thermodynamically  favorable 
in  acidic  media  at  electrode  potentials  below  0.00  V  vs.  RHE  [33]. 

H2^2H++2e~  £$6  =  0.00  V  vs.  RHE  (R6) 

Several  observations  support  the  conclusion  that  reaction  (R6)  is 
contributing  to  the  cathode  current  density.  First,  the  cathode  po¬ 
tential  falls  below  the  limit  of  water  stability  in  acidic  media  (line  a\ 
in  Fig.  4  [34])  at  high  current  density,  indicating  reaction  (R6)  is 


Fig.  4.  Anode  and  cathode  potentials  measured  vs.  Ag/AgCl  reference  electrodes  and 
then  corrected  to  RHE. 
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thermodynamically  favorable.  Second,  an  ex  situ  experiment  in 
which  the  cathode  is  placed  in  a  bath  of  oxidizer  solution  and 
cycled  through  the  same  potential  range  as  in  the  polarization 
curve  measurements  shows  increased  current  density  at  potentials 
below  0.0  V  vs.  RHE  and  rapid  gas  production  (presumably  H2)  at 
potentials  near  -0.4  V  vs.  RHE.  A  plot  of  current  density  measured 
in  the  ex  situ  experiment  is  included  as  Supplemental  information. 
Finally,  extending  the  20  mM  BH4  polarization  curve  to  0.3  V  shows 
an  increase  in  current  density  (Fig.  5),  commensurate  with  onset  of 
reaction  (R6)  in  oxidizer  solution  with  ample  H+.  Reaction  (R6)  is 
included  in  the  calibrated  model  because  of  these  observations. 

Onset  of  the  H202  transport  limit  and  reaction  (R6)  can  be 
observed  separately,  in  contrast  to  the  10  and  20  mM  BH4  curves 
of  Fig.  3,  by  measuring  a  polarization  curve  with  a  membrane 
having  lower  ohmic  resistance.  The  lower  membrane  resistance 
does  not  affect  the  H202  transport  limit,  but  does  shift  the  onset  of 
reaction  (R6)  to  a  lower  cell  potential.  Since  membrane  ohmic 
resistance  increases  with  time  in  DBFC  experiments  [31,35  ,  this 
can  be  accomplished  by  comparing  measurements  with  a  fresh 
membrane  to  measurements  with  the  original  membrane  that  was 
used  in  prior  experiments.  A  fresh  vs.  original  membrane  com¬ 
parison  is  shown  in  Fig.  5,  where  the  decrease  in  cell  potential  at 
-21  mA  cm-2  for  the  fresh  membrane  is  due  to  the  H202  transport 
limit. 

Measured  electrode  potentials,  shown  in  Fig.  4,  can  explain  the 
difference  between  the  theoretical  3.01  V  OCP  for  reaction  (Rl)  and 
the  observed  -1.64  V  OCPs  in  Fig.  3.  The  measured  anode  OCP  is 
1.24  V  vs.  RHE,  which  is  approximately  equal  to  the  expected  value 
E^2.  At  the  cathode,  however,  the  measured  OCP  is  0.4  V  vs.  RHE 
which  is  much  lower  than  the  expected  value  E^3.  The  departure 
from  theoretical  OCP  therefore  occurs  at  the  cathode,  and  is  likely 
due  to  rapid  heterogeneous  H202  decomposition  via  reaction  (R5). 
Consumption  by  reaction  (R5)  depresses  the  H202  concentration  at 
the  cathode-channel  interface,  and  according  to  the  Nernst  equa¬ 
tion,  makes  the  equilibrium  potential  for  reaction  (R3)  less  positive. 
Rapid  gas  production  (presumably  02)  at  OCP  is  confirmed  by 
placing  the  cathode  in  a  bath  of  oxidizer  solution. 

4.2.  Model  calibration 

The  anode  reaction  mechanism  consists  of  reactions  (R2)  and 
(R4),  which  have  rates  given  by  Eqs.  (7)  and  (8).  In  Eq.  (7),  the 
first  order  dependence  with  respect  to  CB h4  is  supported  by 
experimental  evidence  [36  ,  and  CBq2  is  included  to  maintain 


thermodynamic  consistency,  despite  the  rate  being  relatively 
insensitive  to  CBq2  .  Likewise,  the  rate  for  reaction  (R4)  is  effectively 


irreversible  [36],  so  Eq.  (8)  lacks  a  reverse  rate. 

rR2  =  ka.R2  Qh4  A(/>a  -  kcR 2 CBq2  e_^,i?2/A0a  ( 7 ) 

rR4  =  kf,R4C bh4  (8) 

The  cathode  reaction  mechanism  consists  of  reactions  (R3),  (R5) 
and  (R6)  which  have  rates  given  by  Eqs.  (9)— (11). 

rR3  =  ka,R 3e^a^^c  -  kc R3 Cj-j2 o2 e^cR3^A^c  (9) 

rR5  =  ^f,R5^H202  (10) 

rR6  =  -  kCtR6CH2e~ (1 1 ) 


The  second-order  dependence  on  CH+  and  ne  =  2  in  Eq.  (11) 
imply  that  the  rate  limiting  step  for  reaction  (R6)  is  H+  accepting 
an  electron,  which  must  occur  twice.  The  concentration  of  H+  was 
included  in  Eq.  (11),  despite  H+  being  in  excess,  because  it  was  the 
only  reactant  in  the  rate  equation  and  omitting  it  would  have 
permitted  a  “runaway”  reaction  with  no  H+  present.  While  this 
cannot  occur  in  the  real  cell,  it  could  have  caused  numerical 
problems  when  solving  the  model. 

The  parameter  fitting  process  yields  the  values  in  Table  1.  The 
predicted  polarization  curves  are  compared  to  the  measured  curves 
in  Fig.  3,  which  shows  good  agreement  for  the  BH4  limited  stoi¬ 
chiometries  (R2  values  of  0.920,  0.986,  and  0.996  for  1,  2.5  and 
5  mM  BH4  curves,  respectively)  and  reasonable,  but  less  accurate, 
predictions  for  the  H202  limited  stoichiometries.  The  fitted  rate 
parameters  for  reactions  (R2)  and  (R3)  are  not  directly  comparable 
to  the  values  reported  by  Santos  [37]  and  Finkelstein  [38,39], 
because  these  authors  reported  rates  in  terms  of  overpotential 
rather  than  the  electrode-interface  potential  differences  used  here. 
Furthermore,  Finkelstien  calculated  overpotentials  by  assuming  the 
equilibrium  potential  for  each  reaction  was  equal  to  the  observed 
onset  potential.  Nevertheless,  the  fitted  rate  parameters  are  similar 
to  the  reported  values.  Finkelstein  reported  /<a^2  =  6.2  x  10~4  m  s^1 
and  f]a,R2  =  0-2  on  Au  for  5  mM  NaBH4  in  1  M  NaOH,  and 
kc,R3  =  8  x  1CT3  m  s_1  and  (3c,r3  =  0.45  on  Pt  for  5  mM  H202  in  0.5  M 
H2S04. 

The  fit  between  the  predicted  and  measured  polarization  curves 
with  20  mM  BH4  (dash  lines  in  Fig.  3)  could  not  be  improved  with 
the  reaction  rate  for  (R6)  written  as  in  Eq.  (11).  The  discrepancy  at 
low  cell  voltage  may  be  due  in  part  to  the  lack  of  competition  for 
catalyst  surface  sites  in  Eq.  (11 ).  The  adsorption  of  OH  plays  a  major 
role  in  H202  reduction  [40]  and  the  onset  of  significant  H+  reduc¬ 
tion  could  be  delayed  to  a  lower  cathode  potential  by  competition 
with  surface  adsorbed  OH  occupying  sites  on  the  Pd:Ir  surface, 
although  this  study  yielded  no  direct  evidence  to  support  this 
hypothesis. 

Table  1 

Fitted  rate  parameters  for  reactions  R2  through  R6. 


Anode  parameters  Cathode  parameters 


1<a,R2 

9.25  x  10  3  m  s  1 

kc,R3 

7.54  x  10"3  m  s  1 

Pcl,R  2 

0.098 

fic,R3 

0.455 

l<fR4 

3.09  x  loams’1 

kf,R5 

6.34  x  10”4  m  s’1 

Za 

2.73 

kf,R6 

1.19  x  10-9  m4  kmol 

fic,R6 

0.141 

4.11 

Fig.  5.  Measured  20  mM  BH4  polarization  curves;  “Original  Membrane”  is  replotted 
here  from  Fig.  3,  but  with  a  point  at  0.3  V  cell  that  was  omitted  in  Fig.  3. 
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Comparing  the  model  predictions  to  measured  transport- 
limited  current  densities  indicate  that  the  model  correctly  pre¬ 
dicts  the  rates  of  reactant  transport  to  the  electrodes.  The  model 
and  experiments  agree  at  low  cell  potential  for  the  for  the  BH4 
limited  stoichiometries  in  Fig.  3  (1,  2.5,  and  5  mM  BH4  curves) 
showing  that  the  model  accurately  predicts  the  rates  of  BH4; 
transport  to  the  anode.  Turning  off  reaction  (R6)  reveals  that  the 
model  predicts  an  H2O2  transport  limited  current  density  of 
~21  mA  cnrT2  (solid  10  and  20  mM  BH4  curves  in  Fig.  3),  in 
agreement  with  the  measurements  shown  in  Fig.  5  (“fresh  mem¬ 
brane”),  showing  the  model  accurately  predicts  the  rate  of  FI2O2 
transport  to  the  cathode. 

4.3.  Results  from  the  calibrated  model 

A  series  of  simulations  with  the  calibrated  model  were  used  to 
examine  the  influences  of  parasitic  reactions,  cell  geometry,  and 
operating  conditions  on  F1202-DBFC  performance.  H2O2-DBFC  per¬ 
formance  in  each  simulation  was  characterized  by  7]ruj0  ricej0  77th. 
and  power  density.  The  baseline  case  includes  50  mM  NaBH4/2  M 
NaOH  fuel,  250  mM  H2O2/I  M  H2SO4  oxidizer,  and  flow  rates  of 
10  mL  min-1.  These  conditions  were  chosen  for  three  reasons.  First, 
the  baseline  concentrations  are  higher  than  in  the  calibration  ex¬ 
periments,  and  therefore  closer  to  the  concentrations  that  might  be 
used  in  a  practical  F1202-DBFC,  but  still  low  enough  for  limited  gas 
production,  which  preserves  the  validity  of  the  incompressible  and 
laminar  flow  simplifications.  Second,  the  baseline  concentrations 
have  nearly  1:1  stoichiometry  and  should  therefore  reveal  losses  at 
both  electrodes.  Finally,  the  baseline  H2O2  concentration  is  suffi¬ 
cient  to  avoid  significant  H+  reduction  at  the  cathode,  so  baseline 
simulations  address  operating  conditions  for  which  the  calibrated 
model  is  most  accurate. 

4.3 A.  Effects  of  parasitic  reactions 

Losses  due  to  parasitic  reactions  are  manifested  in  two  ways: 
loss  of  reactants  to  gas  production  (lower  77 ce,k),  and  lower  cell 
potential  via  the  Nernst  equation  as  reactant  concentrations  near 
the  electrodes  are  depressed  (lower  77th).  Cell  potentials  predicted 
by  the  calibrated  model  with  BH4  hydrolysis  (reaction  (R5))  and 
H2O2  decomposition  (reaction  (R6))  turned  on  are  lower  than  those 
predicted  with  these  reactions  turned  off,  as  shown  in  Fig.  6.  Both 
curves  in  Fig.  6  are  plotted  at  100  mA  cm-2  current  density  to  make 
ohmic  losses  comparable,  evidenced  by  both  curves  having  the 
same  electric  potential  gradient  in  the  membrane.  The  predicted 
power  density  with  reactions  (R5)  and  (R6)  is  80  mW  cm-2, 


Fig.  6.  Electric  potential  profiles  across  the  cell  predicted  by  the  calibrated  model,  at 
the  midpoint  (25  mm  from  the  inlets),  for  the  baseline  case  at  100  mA  cnrr2  current 
density.  The  two  curves  compare  results  with  reactions  R4  and  R5  turned  on  and  off. 


whereas  the  predicted  power  density  without  them  is  nearly  twice 
as  high  at  152  mW  cm-2. 

The  predicted  effects  of  gas  production  and  depressed  cell  po¬ 
tential  on  the  anode  are  plotted  in  Fig.  7,  which  shows  that  the  fuel 
conversion  efficiency  77C)BH-  has  a  peak  due  to  opposing  trends  in 
?7ce,BH4  and  77th-  The  cell  power  density  is  included  in  Fig.  to  show 
that  peak  fuel  conversion  efficiency  coincides  with  peak  power.  The 
most  energy  can  therefore  be  extracted  from  a  limited  fuel  supply 
by  operating  the  cell  at  peak  power,  in  contrast  to  most  other  fuel 
cells,  such  as  PEMFCs,  which  do  not  suffer  from  parasitic  side  re¬ 
actions.  The  most  efficient  operating  point  for  a  PEMFC  is  near  OCP 
because  thermodynamic  efficiency  is  greatest  near  OCP.  The  model 
predicts  that  peak  oxidizer  conversion  efficiency  also  coincides 
with  peak  power  density. 

Predicted  rice  increases  with  current  density,  and  two  explana¬ 
tions  are  plausible.  First,  increasing  current  density  leads  to  lower 
cell  potential  and  therefore  smaller  activation  energy  barriers  for 
charge  transfer  reactions.  This  favors  charge  transfer  reactions, 
which  then  outcompete  chemical  parasitic  reactions  for  reactants 
at  the  electrode-channel  interfaces.  Second,  high  current  density 
leads  to  lower  reactant  concentrations  at  the  electrode-channel 
interfaces  as  stronger  concentration  boundary  layers  develop  in 
response  to  higher  reactant  consumption  rates.  The  rates  of  para¬ 
sitic  chemical  reactions  are  more  sensitive  to  reactant  concentra¬ 
tion,  and  thus,  they  slow  more  than  the  charge  transfer  reactions  in 
response  to  the  reduced  concentrations. 

4.3.2.  Influence  of  cell  geometry 

The  rates  of  gas  production  by  the  parasitic  reactions  (R5)  and 
(R6),  and  the  net  current  density,  are  plotted  with  respect  to  dis¬ 
tance  from  the  inlets  in  Fig.  8(a).  Gas  production  rates  and  current 
density  both  decrease  with  distance  from  the  inlet,  as  concentra¬ 
tion  boundary  layers  develop  and  limit  concentrations  at  the 
electrode-solution  interface.  The  predicted  gas  production  rates  are 
credible  because  they  depend  on  fitted  reaction  rate  parameters 
and  aqueous  species  transport  for  which  the  diffusion,  migration, 
and  advection  fluxes  can  be  estimated  accurately.  The  predicted 
rates  of  H2  and  O2  transport  in  the  channel,  however,  are  likely  in 
error  because  the  predicted  concentrations  exceed  the  saturation 
limits  for  these  species  in  water.  The  saturation  limit  of  H2  in  pure 
water  at  298  K  is  7.8  x  10-4  M,  and  the  saturation  limit  of  O2  in  1  M 
H2S04  is  9.3  x  10-4  M  [41  .  Higher  concentrations  imply  bubble 
formation  and  because  the  model  does  not  include  multiphase 


Fig.  7.  Effective  fuel  conversion  efficiency,  plotted  with  coulombic  and  thermodynamic 
efficiencies  to  show  how  conversion  efficiency  is  the  result  of  opposing  trends,  and 
plotted  with  power  density  to  show  that  peak  conversion  efficiency  and  peak  power 
density  coincide. 
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(a) 


Fig.  8.  (a)  Predicted  current  density  and  H2  and  02  production  rates  with  respect  to 
position  in  the  channel,  (b)  Influence  of  the  H2  and  02  production  rates  on  the 
coulombic  efficiency.  Both  plots  show  the  50  mM  NaBH4/250  mM  H202  baseline  case  at 
1.1  V  cell  and  10  mL  min-1. 

flow,  predicted  gas  concentrations  are  useful  only  as  indicators  of 
likely  bubble  formation.  Most  bubbles  should  form  near  the  inlets 
due  to  the  preponderance  of  gas  production  there;  this  insight 
could  be  useful  for  the  design  of  graded  catalyst  structures  which 
discourage  bubble  adhesion  near  the  inlets  and  emphasize  reactant 
access  to  the  catalyst  at  points  further  from  the  inlet. 

Fig.  8(b)  plots  predicted  rice  with  respect  distance  from  the  inlets 
and  shows  that  for  these  operating  conditions,  the  cathode  delivers 
electrons  to  H202  more  effectively  than  the  anode  extracts  elec¬ 
trons  from  BH4.  More  than  50%  of  the  H202  flux  to  the  cathode  is 
electrochemically  reduced  but  less  than  40%  of  the  BH4  flux  to  the 
anode  is  electrochemically  oxidized.  Low  overall  ^ce  suggests  that 
future  research  should  aim  to  shift  the  relative  rates  of  reaction 
more  in  favor  of  charge  transfer  reactions.  One  approach  is  to 
develop  novel  catalyst  materials  and  morphologies  [19],  but 
another  is  to  choose  operating  conditions  which  favor  charge 
transfer  reactions.  Among  the  operating  conditions  that  can  be 
changed  to  favor  charge  transfer  reactions  are  the  cell  voltage,  inlet 
concentrations,  and  fuel  and  oxidizer  flow  rates. 

4.3.3.  Influence  of  cell  voltage 

Cell  voltage  can  be  controlled  to  achieve  higher  ^ce.  The 
coulombic  efficiencies  of  the  anode  and  cathode  for  the  baseline 
case  are  plotted  with  respect  to  cell  potential  in  Fig.  9,  which  shows 
that  lower  cell  potentials  increase  rice  at  both  electrodes.  In  effect, 
lower  cell  potential  decreases  the  activation  energy  barriers  to 


Fig.  9.  Predicted  power  density  and  coulombic  efficiency  for  50  mM  BH4/250  mM 
H202  with  10  mL  min-1  flow  rates. 


charge  transfer  reactions  and  therefore  increases  their  rates.  The 
rates  of  parasitic  side  reactions  do  not  depend  on  potential,  and  so 
their  rates  decline  relative  to  charge  transfer  reactions  for  reactants 
at  the  electrode  interfaces. 

Power  density  is  included  in  Fig.  9  to  show  the  relationship 
between  ^ce  and  the  desirable  operating  envelope  of  the  fuel  cell. 
Since  VCeii  falls  with  increasing  power  density,  ^ce  rises  with 
increasing  power  density.  Gains  in  rice  beyond  the  peak  power  point 
are  undesirable,  however,  because  they  coincide  with  falling  7]th  and 
are  accompanied  by  greater  total  reactant  fluxes  to  the  electrodes. 
Beyond  the  peak  power  point,  any  gains  in  ^ce  are  offset  by  overall 
efficiency  losses  and  greater  reactant  consumption  rates,  thus  the 
maximum  desirable  ^ce  is  at  peak  power. 

Coulombic  utilization,  riCUl  is  greater  at  the  anode  because  the 
diffusivity  of  BH4  (2.42  x  10-9  m2s-1  [42])  is  larger  than  H202 
(1.19  x  1CT9  m2s^1  [24  )  and  because  BH4  transport  is  aided  by 
migration  [21  .  The  fuel  and  oxidizer  coulombic  utilizations  both 
show  trends  with  respect  to  Vceii  which  are  similar  to  the  trends  in 
Fig.  9  for  rice.  The  generally  slow  reactant  transport  through 
aqueous  solution  leads  to  low  ricu  (<4%)  for  both  sides  of  the  cell 
under  these  operating  conditions. 

4.3.4.  Influences  of  fuel  flow  rate  and  inlet  concentration 

H202-DBFC  analysis  without  parasitic  reactions  has  shown  that 
raising  reactant  flow  rates  improved  power  density  and  ?7th  through 
increased  convection  mass  transport  in  the  channels  [21  .  The  same 
is  true  when  parasitic  reactions  are  included,  but  with  a  caveat.  The 
higher  rates  of  transport  to  the  electrodes  also  lead  to  higher  rates 
of  gas  production.  The  influence  of  fuel  flow  rate  on  concentration 
boundary  layer  development  is  shown  in  Fig.  10,  where  the  BH4 
concentration  boundary  layer  at  10  mL  min-1  fuel  flow  rate  is 
compact  and  the  boundary  layer  at  1  mL  min'1  is  beginning  to 
envelop  the  entire  channel.  The  compact  boundary  layer  is  the 
result  of  convection  augmenting  the  rate  of  diffusion  mass  trans¬ 
port  to  the  anode. 

Variation  in  peak  power  density  among  curves  in  Fig.  11(a) 
shows  how  greater  power  density  incurs  a  lower  upper  bound  on 
coulombic  efficiency.  As  fuel  flows  more  quickly,  BH4:  concentration 
at  the  anode  interface  rises  in  response  to  more  facile  transport 
from  the  bulk.  Higher  BH4  concentration  favors  hydrolysis,  so 
coulombic  efficiency  falls.  Peak  power  density  approaches  the 
horizontal  axis  with  decreasing  flow  rate  until  it  intercepts  the  axis 
at  zero  flow  rate.  The  maximum  theoretical  coulombic  efficiency  for 
these  operating  conditions  is  -73%  and  occurs  in  the  zero  flow  rate 
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Fig.  10.  Predicted  BH4  concentration  in  the  fuel  channel  for  the  baseline  case  at  1.1  V 
for  (a)  10  mL  min-1  and  (b)  1  mL  min-1  flow  rates.  Both  plots  share  the  same  color 
map.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 


case,  where  diffusion  and  migration  alone  dictate  the  rates  of 
transport  to  the  anode.  The  power  density  shows  a  different  rela¬ 
tionship  with  respect  to  coulombic  utilization.  In  Fig.  11(b),  peak 
power  density  asymptotically  approaches  zero  as  the  coulombic 
utilization  approaches  100%. 


Fig.  11.  Predicted  (a)  coulombic  efficiency  (b)  and  coulombic  utilization.  Both  plots 
share  fuel  flow  rates  from  1  to  15  mL  min-1,  oxidizer  flow  rate  of  10  mL  min  \  and 
50  mM  BH4/25O  mM  H202. 


The  influence  of  inlet  concentration  on  H2O2-DBFC  performance 
is  examined  with  respect  to  the  concentration  of  BF14 .  Changing  the 
inlet  concentration  revealed  a  relationship  between  power  density 
and  coulombic  efficiency  which  mirrors  that  of  the  fuel  flow  rate. 
Increases  in  fuel  flow  rate  and  inlet  concentration  both  shift  the 
relative  rates  of  reaction  at  the  anode  in  favor  of  H2  production  by 
increasing  BH4  concentration  at  the  anode  interface.  The  coulombic 
utilization  trend  with  respect  to  BH4  concentration  differs  from  the 
trend  with  respect  to  flow  rate.  Slower  flow  rates  increase  residence 
times  in  the  F1202-DBFC  while  diffusion  and  migration  continue  to 
transport  BH4  to  the  anode,  leading  to  the  asymptotic  approach  to 
zero  power  density  at  100%  coulombic  utilization  seen  in  Fig.  11(b). 
Lower  inlet  concentrations  decrease  the  rate  of  transport  by  all 
three  processes  (migration,  diffusion,  and  convection),  leading  to  a 
theoretical  upper  bound  of  -4.75%  coulombic  utilization  at  zero 
inlet  concentration. 


5.  Conclusions 

We  show  that  global  reactions  and  rates,  if  properly  calibrated, 
can  describe  borohydride  oxidation  and  hydrogen  peroxide 
reduction  with  sufficient  accuracy  for  H2O2-DBFC  design  analysis. 
The  finite  volume  approach  applied  to  modeling  transport  in  a 
H2O2-DBFC  accurately  predicts  the  fluxes  of  ions  and  neutral  spe¬ 
cies  in  the  reactant  solutions,  without  the  complexity  of  a  non-ideal 
electrolyte  model,  for  the  operating  conditions  examined  in  the 
experiments. 

A  numerical  model  calibrated  to  a  simple  H202-DBFC  shows  that 
the  limited  diffusivities  and  mobilities  of  reactants  in  aqueous  so¬ 
lution  are  responsible  for  the  development  of  compact  boundary 
layers  and  low  single-pass  reactant  utilization.  The  relative  rates  of 
charge  transfer  reactions  and  parasitic  reactions  are  shifted  in  favor 
of  charge  transfer  by  lower  reactant  concentrations  near  the  elec¬ 
trodes  and  larger  activation  overpotentials.  These  conditions  can  be 
achieved  by  lower  reactant  inlet  concentrations,  slower  bulk  flows 
in  the  channels,  and  operating  the  cell  at  high  current  density  (low 
cell  voltage).  The  opposite  trends  in  cell  coulombic  efficiency  and 
thermodynamic  efficiency  with  respect  to  current  density  lead  to  a 
peak  overall  efficiency  point  which  (for  the  cell  geometry  and 
operating  conditions  examined  in  this  study)  coincides  with  peak 
power  density.  The  coincidence  of  peak  conversion  efficiency  and 
peak  power  distinguishes  H2O2-DBFCS  from  fuel  cell  technologies 
that  are  less  strongly  affected  by  parasitic  electrode  reactions.  The 
insight  into  H2O2-DBFC  behavior  provided  by  the  numerical  model 
illustrates  the  utility  of  such  models  for  cell  design  and  selection  of 
appropriate  operating  parameters  when  complex  chemistry  and 
transport  processes  complicate  performance  prediction. 

In  future  work,  the  demonstrated  modeling  approach  could  be 
applied  to  other  electrochemical  systems,  including  alternative 
H2O2-DBFCS  utilizing  alkaline  oxidizer  solutions  and  anion- 
exchange  membranes.  The  transport  and  membrane  sub-models 
would  be  modified  to  include  alternative  species  (HO2  in  the  case 
of  an  alkaline  oxidizer  solution  and  OH-  in  the  case  of  an  anion- 
exchange  membrane)  and  the  cathode  reaction  rate  sub-model 
would  be  adapted  to  a  different  reaction  mechanism.  The  govern¬ 
ing  equations,  solution  approach,  and  calibration  methods  would 
remain  the  same. 
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Appendix  A.  Nomenclature 

a  anodic  parameter  (subscript) 

j3  electrochemical  reaction  symmetry  factor 

c  cathodic  parameter  (subscript) 

Q  concentration  of  species  k 

standard  half-cell  potential  for  reaction  R 
F  Faraday's  constant 

f  FIRigT 

i  current  density 

J  mass  flux 

J  k,e lec  flux  °f  species  k  to  an  electrode 
J  k  mole  flux  of  species  k 

kq  rate  constant  for  reaction  q 

£  roughness  factor 

LCeii  length  of  flow  channels  from  inlet  to  outlet 
7]rUik  reactant  utilization  with  respect  to  species  k 
Tlth  thermodynamic  efficiency 

Vcujc  coulombic  utilization 

rice  coulombic  efficiency 

rjCtk  conversion  efficiency  with  respect  to  species  k 

ne  number  of  e_  in  rate  limiting  reaction  step 

neiec  number  of  e_  transferred  per  mole  of  reactant  consumed 

</>  electric  potential 

P  pressure 

p  solution  mass  density 

pc  net  charge  density 

Rig  ideal  gas  constant 

rq  rate  of  reaction  q 

T  temperature 

uk  mobility  of  species  k 

vx,vy  velocity  in  x-  and  y-  directions 

vk,q  stoichiometric  coefficient  of  species  k  in  reaction  q 

V  volumetric  flow  rate 

x  coordinate  axis  parallel  to  membrane 

y  coordinate  axis  perpendicular  to  membrane 

Yk  mass  fraction  of  species  k 

zk  charge  number  of  species  k 

jk  activity  coefficient  for  species  k 

Appendix  B.  Supplementary  data 

Supplementary  data  related  to  this  article  can  be  found  at  http:// 
dx.doi.org/10.1016/j.jpowsour.2014.07.139. 
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